Maximal superintegrability of Benenti systems 
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Abstract 

For a class of Hamiltonian systems naturally arising in the modern theory of separation of 
variables, we establish their maximal superintegrability by explicitly constructing the additional 
integrals of motion. 

For a Hamiltonian dynamical system on a 2n- dimensional phase space to be completely integrable, 
the classical Liouville theorem requires n involutive integrals of motion. Solving the equations of motion 
then amounts to n quadratures, but performing the latter still can be a challenge. However, if there 
are m > n integrals of motion, only 2n — m quadratures are required. In particular, for the so-called 
maximally superintegrable systems (i.e., those possessing 2n — 1 integrals of motion) we need just one 
quadrature. 

Superintegrable systems on two- and three-dimensional Euclidean spaces and spaces of constant cur- 
vature were extensively studied in the literature, see e.g. CP-[ZI and also ^ for the case of nonconstant 
curvature. On the other hand, there is not much known about superintegrability in higher dimensions. 
Remarkable exceptions include n-dimensional Winternitz-Smorodinsky model PP and its generalizations 
P, Calogero-Moser-Sutherland system jTUllII] and the systems with isochronic potentials [12] and mod- 
ified Coulomb potential Notice that for the majority of these systems the separation of variables 
takes place in the original coordinates and the Hamiltonians themselves are additively separable. 

In this letter we consider the Stackel systems from the class described by Benenti |14[ [T3] , so the 
corresponding Hamiltonians are quadratic in momenta. We show that some of these systems (namely 
those associated with flat or constant curvature metrics) are maximally superintegrable for arbitrary 
dimension n of the configuration space, and it seems plausible that they are multiseparable too. 



Consider an n-dimensional manifold M with local coordinates q^, 
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If we interpret G as a contravariant metric on M, then L is a special conformal Killing tensor for G in 
the sense of Crampin and Sarlet |TE] (and moreover, this holds for any metric of the form UG, -i G Z 
[T7]). and hence for any integer i the quantities 
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are in involution with respect to the canonical Poisson bracket in (p, g)-variables: {q^,Pj} = Sj, i,j = 
1, . . . ,n. Here (g^, . . . , g", pi, . . . , p„) are local coordinates on the cotangent bundle T*M and the Killing 
tensors Kr are constructed from L as follows 

r 

K, = I, ir,+i = J]g'=L'-^ r = l,...,n-l, 

where we set for convenience = 1 and I stands for the nxn unit matrix. Notice that are coefficients 
of the characteristic polynomial of L: 

det(e/ -L)= C + q'C + - + g"- (2) 

The Hamiltonians Ei r belong to the family of Stackel separable systems, and the separation of 
variables is achieved by passing to the coordinates (/i, A) related to {p, q) by the formulas 

n 

P. = - 5^(A^)"->fc/Afc, q' = (-1)V,(A), z = l,...,n, 

k=l 

where is an i^^ symmetric polynomial in the variables A^, . . . , A" (uo = 1, o"i = ^17=1 ■ ■ ■ I'^n = 
AiA2...A")andA, = n;=i,y.(A^-A^). 

The main advantage of using the 'nonstandard' variables (p, q) is that Ei^^ are polynomial in g's for 
i > 0, so it is natural to search for additional integrals of motion that also are polynomial in p's and g's. 

For 2 = 0, . . . , n the (contravariant) metrics {GiY^ = [UCy^ are flat and we have 

n— t— 1 71— i i k 

EiA = 2 X] 9''^ X] PjPn-i+k-j+1 ~ 2 ^ S Pn-i+jPn-i+k-j+l- 

k=0 j=k+l k=l j=l 

The metric Gn+i has nonzero constant curvature and we obtain 

n 

En+1,1 = I E {(I'Q^ - (f^^)PiPj. 

where we assume that = for > n. 
We readily flnd that 

{Ei,i,p„„i} = 0, 2 = 0, ...,n-l, 
{Ei,i,g>„_i+2} = 0, z = 2, + 

{Ei,i,g^} = -pn-u i = 0,...,ra-l. 

Thus g"~* is a cyclic coordinate for Ei^i, i = 0, . . . ,n — 1, and Fi^r = {Ei^r,Pn-i} for r = 2,...,n 
provide {n — 1) additional, quadratic in momenta, integrals for the dynamical systems associated with 
the Hamiltonians Ei^i, i = 0, . . . ,n — 1. 
Consider s functions of the form 

n 

U= Y,Al'{q)p,p, + aiq). 
j,k=i 

It is easily seen that if the matrices Ai are linearly independent over the fleld of all locally analytic 
functions of g^, . . . , g",pi, . . . ,Pn, then fi are functionally independent. 

Using this result we readily flnd that (2n — 1) quantities Ei^s, s = 1, . . . , n, and Fi^r, r = 2, . . . , n 
are functionally independent for any i = 0, . . . ,n — 1, and hence the Hamiltonians Ei^i are maximally 
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superintegrable for all these i. The same is true for i = n and i = n + 1, but in this case the additional 
integrals have the form Fi^r = {Ei^r, (fPn-i+2}, r = 2, . . . ,n, and they again are quadratic in momenta. 

This picture can be extended to the case of nonzero potentials. Namely, consider the basic separable 
potentials given by the recursion relations IT^ 

v;("+i) = v;(-) + \//iV/™\ m = 1, 2, . . . , Vr^'^ = -g^ 

Vr^'^ = 0, (3) 

K^-™-!) = v}:r^ + \/^'Vi-'"\ m = 1, 2, . . . , K^-') = -f-^/q^. 

As the potentials v}^"^ are independent of for m = j — n, . . . , j — 1, the above analysis can be 
generalized to yield the following result: 

Theorem 1 For any natural n > 2 the dynamical systems with the Hamiltonians Hj^f = Ei^i + v}'''^ are 
maximally superintegrable for alii = 0, ... , n—1, k = —i, n—l—i and i = n, n+1, k = 2—i, ...,n+l—i. 

(k) (k) 

Namely, besides n involutive integrals HlJ = Ei^r + K , t = 1, . . . , n, the dynamical system associated 
with h[^^ has {n — 1) additional integrals F-'l^, s = 2, . . . ,n, of the following form: 



«^ ^ff = {Hl^J,Pn-i}, s = 2, . . . ,n, for i = 0, . . . ,n - 1, k = -i,...,n- 1-i, 
b) F^'^s = {h[% q"Pn-i+2} fori = n,n + l, k = 2 - i, ...,n + 1 - i, 
and these {2n — 1) integrals (H^^J , r = 1, . . . ,n, and F^^] , s = 2, . . . ,n) are functionally independent. 
Let us illustrate this theorem by a nontrivial example. First, let n = 4, i = A, k = —2. The functions 

Tr( — 2) -^122 4 -^'^24 i'^// 4\2 

^4,1 = Pi -Q P1P2 - q PiPs - q piPa - ^q P2 - g P2P3 + q /{q > , 

-P2qYp3 - \qW + \{q')W + q'q%P^ + \{q')W + q'qV{q'f - l/q\ 
11 1 

h'i^^ = --q^pi' - g^piP2 + -q^q^P2^ + {q^fp2P2, + q^q^P2PA - -q^q'^P2^ + q^q^p^'^ 

Hqypm-q'/q' + qV/iqr, 
Hif = -\qW + \q'qW + qWp2P, + {q'fp2P, + ^(g^)V3^ " ^V?' + {q^Yliq'? 

are in involution by construction. Moreover, h\ commutes with q'^p2 and thus by Theorem 1 the 



^^4,2 = ~2^^^^^ ~ ^^^1^2 - q^piPz + i^{q^fp2^ + q'^q^P2P3 + q^q^P2PA - ■^q^q^P2 



quantities 

7("2) — -^^4 2 p(-2) _ 1 4 3 2 I ^„4\2^ ^ ^3 / 4 p(-2) _ ^ (^i\2 2 , 

4,2 - 2^ -^4,3 - i^q q P2 + \q ) P2P3~ q /q , ^4,4 -2^q>P2 + ^^ 



are additional integrals of motion. As ifj ij|2^^ -^43^''? -^4 4^''5 -^4 2^'*5 -^4 3^^ -^44^"* functionally 
independent, the dynamical system associated with ^ is maximally superintegrable. Notice that all 
(2n — 1) integrals in this case are quadratic in momenta. 

For all natural n > 2 we have dv}"^^ /dq^ = — 1 for m = j and dv}"^^ /dq^ = 2q^ for m = j + 1, 
whence {hI}~'\ - gi} = and {if^""*^^^ ^Pl-i + (q^?} = for z = 0, . . . , n - 1. Therefore the 
dynamical systems associated with H^^~^^ and ifj^""*"*"^-* for i = 0, . . . , n — 1 possess (n — 1) additional 
integrals, namely F^"^) = {H^^ \pU - q'}. s = 2,...,n, and F^'-"'^ = {H^'-''\ \pU + {q'f}. 



s = 2, . . . , n, respectively. However, in these cases the additional integrals are cubic in momenta, so we 

(k) (k) 

were unable to prove the functional independence of the sets [H-^ , r = 1, . . . , n, and F^/, s = 2, . . . , n) 
for i = 0, . . . , n — 1 and k = n — i and k = n — i + 1 in full generality so far. Nevertheless we are 
certain that these sets indeed are functionally independent for all these values of i and k, and hence 
the dynamical systems associated with if-*"" *■* and H^^^ are maximally superintegrable. 
For instance, let n = 3, i = 0, k = 4. The functions 

= P^P^ + \pI + ^iP2P3 + \pW + 2g'g' - ^q\q'f + {q^f + {qi)\ 

rrf4) 9 /1\2 "^32 "^122 /1\2S r\2S 

^0,2 = + qipm + giP2 + \q ) P2P3 - ^q P3 + ^q q P3- [q ) q +^q q 

-2q\qY + {q')V, 

-^2 1 2 '^/1\22 12 '^1'^2 

^0,3 = 2^1 + 9 P1P2 + q P1P3 - q P2P3 + ^{q ) P2 + q q P2P3 - ^q q Ps 

+ liqrpl + qH-2qW + q' + {qy) 

again are in involution by construction; Hq'^I also commutes with |p| + (q^)"^, and hence by the above 
so do 

Fo^J = -2gV2-3(g^)V-^pi + 2gV3, 

^0,3 = -q^Pi - 2(g') V - 4q\% - P2pI - ^gVs + 2g'p3 + {q^fp3- 

As Hq'^I, Fq^2 J ^03 functionally independent, the dynamical system associated with Hq^I 

is maximally superintegrable as well. In contrast with the previous example, the additional integrals 
5 -foils' cubic in momenta. 
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